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Abstract. We establish the dynamical attractor behavior in scalar-tensor 
theories of dark energy, providing a powerful framework to analyze classes of 
theories, predicting common evolutionary characteristics that can be compared 
against cosmological constraints. In the Jordan frame the theories are viewed as 
a coupling between a scalar field, and the Ricci scalar, R, F(&)R. The Jordan 
frame evolution is described in terms of dynamical variables m = din F/dln <I> 
and r = — QF/f, where F(&) = df(<fr)/d<&. The evolution can be alternatively 
viewed in the Einstein frame as a general coupling between scalar dark energy and 
matter, j3. We present a complete, consistent picture of evolution in the Einstein 
and Jordan frames and consider the conditions on the form of the coupling 
F and fi required to give the observed cold dark matter (CDM) dominated 
era that transitions into a late time accelerative phase, including transitory 
accelerative eras that have not previously been investigated. We find five classes 
of evolutionary behavior of which four are qualitatively similar to those for f(R) 
theories (which have /3 = 1/2). The fifth class exists only for \/3\ < %/3/4, i.e. not 
for f(R) theories. In models giving transitory late time acceleration, we find a 
viable accelerative region of the (r, m) plane accessible to scalar-tensor theories 
with any coupling, (3 (at least in the range |/3| < 1/2, which we study in detail), 
and an additional region open only to theories with \/3\ < V3/4. 
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1. Introduction 

The current accelerated expansion of the universe has led to a variety of new theories 
in cosmology under a common descriptor of dark energy. The acceleration was first 
observed using data from Supernovae la [TJ |5J |3J HI [5] and has been confirmed by 
complementary observations including WMAP data on Cosmic Microwave Background 
anisotropies (CMB) El[7l[8][9], data from the SDSS on Large Scale Structure formation 
(LSS) [TO! [TT1H2], baryon acoustic oscillations PSES] and weak lensing [T6 | ITT]. 

The simplest way of introducing an acceleration into Einstein's equations is to 
add a cosmological constant so that the decelerating expansion (which we naturally 
expect because of the attractive force of gravity) gets nullified and over-compensated 
for. The cosmological constant though suffers from the problem of fine tuning since its 
energy density is 10 120 times smaller than what we might expect if it related to Planck 
scale phenomena. Other models of dark energy include introducing a scalar field (such 
as minimally coupled quintessence models PHI HS1 1201 HH 122] > coupled quintessence 
models, [23J OH [231 [2BJ [27], K-essence [Ml Hi 132], Chameleon fields [32 [MIES] etc), 
Phantom dark energy [321 GE3 GES] , Chaplygin gases [33 [3S 133 SB S3] , topological 
defects |42j . and many others. A recent review on dark energy can be found in Ref. 

An alternate approach to explaining the observed acceleration, without recourse 
to a new form of energy, is to consider large scale modifications to gravity coming from 
scalar-tensor gravity [4^[11[13[13[13[1E1[10[S3 higher 
order gravity theories E3 [SHI EU E3 E31 EH EH EE] such as modified Gauss-Bonnet 
gravity. One of the earliest modifications [44] added a l/R term to the gravitational 
action. Such modifications to gravity seem to conflict with local solar system tests 
of gravity, however general f(R) gravity models remain interesting because there are 
ways of evading local gravity constraints [3U [321 1331 EZ] • 

Dynamical attractors in dark energy theories provide a powerful way to analyze 
classes of theories. They predict common evolutionary characteristics, with minimal 
dependence on fine-tuning of initial conditions [68], [23] EH [27], with which we can 
critique theories in light of cosmological constraints on the background evolution 
|53l 169] and also the growth of perturbations and implications for large scale structure 

EH- 

In this paper we apply dynamical attractor analysis to general scalar-tensor 
theories. In the Einstein frame theories are described by a coupling between scalar 
dark energy and matter, parameterized by /3, and in the Jordan frame by a general 
coupling between a scalar field and the Ricci scalar, F(Q)R; (3 — and F(Q) = 1 
for General Relativity. We have structured our paper in the following way. Section [2] 



contains the equations for background evolution in the Einstein frame (Sec. 2.1l and 



Jordan frame (Sec. |2.2| for generally coupled models. We discuss the two alternative 



perspectives given by the two frames (Sec. 2.3 1 and the form of the Jordan frame 
potential we consider (Sec. 2.4). In section |3j we set up all of the equations that we 



need to study the critical points of dynamical evolution in the Jordan frame (Sec. 3.1l 



and Einstein frame (Sec. 3.2 1, discuss the mapping of critical points from one frame 
to the other (Sec. |3.3| and how stability is assessed (Sec. 3.4 1. In section [4] we present 



the dynamical attractor solutions (Sec. 4.1l and compare our solutions with those 



previously obtained in the literature for f(R) and non-minimally coupled quintessence 



theories (Sec. 4.2 1. We perform a detailed analysis of the stability of all critical points 



(Sec. 4.3 1 and classify scalar-tensor theories into different classes (Sec. 4.4). We then 
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discuss behavior if we relax the condition of a stable accelerated universe to include 
saddle critical points (Sec. 4.5). Finally we summarize and conclude in section [5] 



2. Cosmic evolution 



2.1. Evolution in the Einstein frame 



We start with the action in the Einstein frame (EF) 

1 ^ 1,^ 
2k 



S = 



d A x 



- L 



d A x 



-2 v /2/3k / 30- 



(1) 



where 



8itG, ^ is a scalar field which is non-minimally coupled to matter through 



a coupling strength (3. We denote all quantities in the EF with a tilde, such as the 
metric, g^u, Ricci scalar, R, scalar potential, V{4>) and the matter Lagrangian, C m . 
The action in the EF is a standard General Relativity action with a non-minimal 
coupling between the scalar field and matter. For a Robertson Walker metric with 
flat geometry described by a cosmic expansion factor, a, evolving in physical time, t, 
the Friedmann, acceleration, scalar field and matter fluid equations are respectively, 

H 2 =^i~p m +t(^\ 2 + v)> (2) 



dt 



dH 
dt 

d^ 
dt 2 

dpm 

dt 



3H^ 
dt 



V lh = 



3-ff (p m + Pm) 




3pm), 



(3) 



(4) 



(5) 



where p m and p m are the EF matter density and pressure, H = dlna/dt and 
V,0 = dV/dcj). In this paper we are interested in the nature of cosmic evolution 
in the matter dominated and dark energy dominated eras. We therefore restrict our 
attention to pressureless matter, with p m = 0. 

Such a parameterization, using p to describe a non-minimal coupling between 
matter and scalar, has been widely discussed in the literature, for example [2"3l |2"4"1 l2"5l 
l26l l53"1 l69l |27j [70j |71]. This can include low temperature neutrinos which behave as 
nonrelativistic matter. The scalar coupling of the form shown in ^ can hence give 
rise to mass varying neutrinos ('MaVaN's) [721 [731 1741 ITT] . 



2.2. Evolution in the Jordan frame 

We can re-express the action (JT|) in the Jordan frame (JF) as a non-minimal 
coupling of the scalar to gravity and minimally coupled matter, through a conformal 
transformation, redefining the metric, 

-?-/273 K /30~ 



g^u = e ~v~>" , " r " r g IM ,- (6) 
This leads to a general Brans-Dicke type action for a scalar-tensor theory with a 
non-canonical kinetic term, 

2 

(V$) 2 - V($) 



S= d 4 i 



-f) 



1 



3(1 - 4/3 2 ) 1 
16k 2 /3 2 F($) 



dF 
d¥ 
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where we have introduced the function of scalar <I>, which is wholly specified by 

4> and /?, 

F($) = e 2 ^^* (8) 

We require .F(<i>) > for the conformal transformation to be well defined. <&(</>) can 
be specified by choosing an explicit form for F(<&) and inverting 

Pressure, density, time, and expansion factor in the two frames are related via 

Pm = F{<S>) 2 p m , Pm = F{$>) 2 p m , dt = F($)- 1 / 2 dt, a = F(^)-^ 2 a, (9) 

with the potential in the JF related to the EF potential by, 

V($) = F(4>) 2 l/(0), (10) 

where again, an explicit form for V(&) can be found given a specific choice of F(<I>) 
and hence <&(</>)• 

The Friedmann, acceleration, scalar field and energy conservation equations for 
the action in the JF are: 



^ (fc + m) _f + <i__|!)|! 

1 2 If 1 HF (1 - 4/3 2 ) F 2 
H ~~2F K Pm ~ 2F + 2^~ 3 16/32 F*' (12) 

F + mF= 4 ^(^p m +W($)-2F^j, (13) 

Pm + 3H Pm - 0. (14) 

where dots denote derivatives with respect to physical time t in the JF. 

2.3. The Einstein frame vs. Jordan frame perspectives 

If a theory can be expressed in a frame in which all matter components are minimally 
coupled, such as in ([T]), then a conformal transformation to the Einstein frame will 
always result in the matter-scalar coupling being the same strength, (3, for all types of 
matter. Some debate persists in the literature, in the case of such identical couplings, 
on whether EF or JF is the 'physical' frame, see for example Refs. [75J |751 [77] 
on this issue. We favor the perspective that neither frame is more physical than 
the other, but that cosmological observations can be viewed more conveniently in 
the JF since interpretation of observations that form the key evidence for dark 
energy, of fluctuations in the CMB, and of redshift measurements in galaxy surveys 
and supernovae, are usually made with the assumption that baryons are minimally 
coupled. When considering the dynamical evolution of scalar-tensor theories therefore 
we consider properties, for example the fractional matter density, equation of state 
and stability of critical points (attractors/repellers), in the JF to be those that are 
most appropriate to compare with observational constraints. 

The Einstein frame perspective is also useful however, in that the evolution 
equations and fixed point analysis are arguably easier analytically in this frame. 
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The Einstein frame description can also be used to encompass theories outside of 
the modified gravity realm. In these theories the scalar is minimally coupled to gravity 
but the coupling strength of the scalar with cold dark matter (/3 C ), baryons ((3b) and 
neutrinos ((3 V ) could be different, for example ( Ell Ell EH1 E31 Ell E01 EI] ) - In 
models in which CDM alone is coupled observational constraints constrain (3 C < 0.1 
[24] . In Chameleon Cosmology models [31] [32] [33] , the scalar fields acquire a mass 
whose magnitude depends on the local matter density and all (3i can be spatially 
varying, and of order unity. 

One of the aims of this paper is to give a coherent picture of dynamical attractor 
analysis in modified gravity theories in both the Einstein and Jordan frames. As will 
come out naturally in our analysis, the choice of frame is really a matter of convenience 
and the physical implications of the analysis, in terms of attractor critical points 
and their stability arises identically out of analyses irrespective of frame choice. We, 
therefore, focus our analysis on scenarios in which the coupling strength is the same 
for all types of matter and in which the Jordan frame quantities are to be compared 
with observational predictions. The Einstein frame analysis however can be equally 
applied to theories in which CDM-scalar couplings are present, and the Einstein frame 
is the physical frame. 



2.4- The form of the scalar potential 

The action in Q allows both a free choice of coupling to gravity F(Q) and scalar 
self interaction potential V(<&). In order to include an important group of extended 
gravity theories, l f(R) theories', and ACDM as classes of models described by our 
analysis, however, we place a restriction on the form of the JF potential. 

f(R) theories, where the Lagrangian contains an arbitrary function of the scalar 
curvature R, can be written as a scalar tensor theory of gravity [75] [50] H3 [ST] , 

S = ~ / d A x^~gf(R) + S m (15) 



2k 2 
1 

2«2 



i 4 Xy 



-g{F(<S>)R - ($F($) - /(<&))] + S rn (16) 



where, 

F($) = d/($)/d$. (17) 

This is equivalent to our general scalar-tensor action ([7]), with (3 — 1/2, i.e. a non- 
dynamical scalar field, and scalar potential, 

1 



V ($) = s-2 [**"(*)-/(*)]• (18) 



By considering actions with potentials of the form ( 18 1, our analysis has f(R) theories 
as its limit when the field $ is non-dynamical. Equally, with this form of the potential, 
general relativity is regained for /($) = $ and (3 = 0. 

In the dynamical analysis described in the following sections, with this restricted 
potential, we find it useful to split the potential in (18) into two components, 

= - i/W' W) = (19) 
and consider their relative importance in the cosmic evolution. This analysis could 
be extended, therefore, to describe evolution of more general actions, without the 
restriction on the potential given in (18}, through writing 



7(*) = _L[*f($) -/($)] + %(*) (20) 
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and considering the relative importance of the extra component, V3(<I>), to the 
dynamical evolution. We leave such an extension to future work. 



3. Dynamical critical points 

3.1. Jordan frame autonomous phase plane equations 

We would like to find the critical points of the dynamical background evolution. This 
equates to finding asymptotic power-law solutions of the form a cx t p , where p is 
a constant and related to the effective equation of state of the system, w e ft, by 
p = 2/3(1 + w efS ). 

Writing the scalar potential as (18 1, and defining the dynamical variables, 



_F___dhi£ _ _J_ _ 
Xl ~ HF~ dlaa' ^ ~ 6FH 2 ' ^"Gif 2 ' [ ' 

the fractional matter density in the JF can be expressed in terms of these dynamical 
parameters through the Friedmann equation (111, as, 

fim = 3 pff2 = 1 ~~ Xl ~ Kx l ~ x 1~ x 3> ( 22 ) 

where K = (1 — 4/3 2 )/16/3 2 . w e ff can also be expressed in terms of the dynamical 
parameters, 

»<.//= -fi^- 1 < 23 > 
_ 4|! + (i^!) T _ + ^ _ ^ Ki; _ (1 _ if) xi _^_ ^ xs (24) 

so that critical points, where w e ff is constant in time, are satisfied by 
0, where x' = dx/dhxa = x/H. 

The closed set of 'autonomous phase plane' equations, x\, x' 2 and x' 3 in terms of 
xi, X2 and 23 are obtained from the JF equations ( 11 ) - ( 14 1, and constraint equation 

x 'i = ~ y [( 3 " ^ ~ ( 3 ~ - ( 3 - ^ 2 )Kx\ + (3 - 12/3 2 )x 2 + (3 - 4/5 2 )x 3 ] 

-4/3 2 (l-a; 1 -ii:.x 2 + 3.T 2 + X3) (25) 
4 = x 2 [(3 + 4/3 2 ) + (2 - 4/3 2 ) xi + (3 - 4/3 2 ) i^ 2 - 3 (l - 4/3 2 ) a? 2 - (3 - 4/3 2 ) 23] 

TO 

4 = 2:3 [(3 + 4/3 2 ) + (1 - Af3 2 ) xi + (3 - 4/3 2 ) ifa: 2 - 3 (l - 4/3 2 ) a* - (3 - 4/3 2 ) 23 
£1X3 



TO 

where m is defined as 



(26) 

"a] 
(27) 



dhiF /din 23 , , 

^[^T-^ + ^ff) ■ (28) 



dln<I> \ din a 

and its explicit evolution depends on the choice of /(<&) 
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3.2. Einstein frame autonomous phase plane equations 



We can define analogous dimcnsionlcss dynamical variables, x, y and z, to those 
defined in Sec. |3.1[ which are simply related to the fractional energy densitites of the 
scalar kinetic and potential energies in the EF, 



where V(4>) 



V6H dt 

VMS 



f 



2k 2 F 2 ' 



3H 2 



v 2 (4>) 



3H 2 



<S>F 
2k 2 F 2 ' 



(30) 



The specific form of the potentials V± and V 2 coming from the conformal 
transformation will depend on the form of /(<&)■ 

We can also define an Einstein frame fractional energy density for matter, and 
using ([2]) write it in terms of the Einstein frame dynamical variables, 



3H 2 



= 1 - (x 2 + y + z) 



(31) 



where < Cl m < 1. 

On differentiating the potentials Vi(<p) and V 2 {4>) with respect to <j> and using the 
scalar field redefinition equation (Tsl) we get, 



2,0 



K(3 



k/3V 2 



2Vi 



m 



1 

m 



(32) 



(33) 



The EF equations ([2| - (|5j, along with the constraint equation (31) can be used 
to write down the system of plane-autonomous equations, 



dx 
din a 

dy 
din a 

dz 
din a 



(3(1 



x 2 



y-z)+2[3(2y + z)--x(l-x 2 + y+~z) 1 (34) 



4f3xz 



8f3xy + 3y(l 



~x 2 



V ~ z), 



1 



A(3xz 1 



35(1 + x 2 



V - z). 



(35) 
(36) 



where all information about the scalar field redefinition </)($) is encoded in the 
parameter m. 

We can find EF evolution where there is power law expansion of the form a oc t p 
where p is constant by solving dx/dlna = dy/dhid — dz/dlnd = 0. As discussed 
previously, if visible matter, and not only CDM, is non- minimally coupled in the 
Einstein frame then Einstein frame is not the one in which cosmological observations 
are usually expressed. However as we discuss in the following section, a mapping exists 
between the autonomous phase planes in each frame that implies that solving one set 
of equations translates into the solutions of the other. 
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3.3. Mapping between the Einstein and Jordan frame phase planes 
There is a direct mapping between the JF and EF variables, 



X\ = 



X2 



A(3x 
l-2(3x 

y 

(1 -2/3i) 2 ' 

z 



1 



.Tl 



y 



2f3 xi -T 

4:X 2 



where 



(1 -2/3£) 2 ' 

din a 1 
dlna ~~ 1 - 2(ix ' 



(Xi-2)2' 

Z Ori-2) 2 ' 



Xl 



(37) 
(38) 
(39) 

(40) 



Writing a point in the JF as x = (x\, x%, x$) and a point in the EF as 
x = (x, y,z), the transformation matrix T defined such that Ax = TAx, where 
Ax = {Ax i, Ax2, Ax^) and Ax = (Ai, Ay, A5), is given by, 



r = 



-4/3 
(1 - 2l3~xY \ 1 i|f 



1 



1-2/35 



1 



r- 1 = 



1 

Sx 2 A 



(xi-2) 






4 



(41) 



(42) 



The transformation from the Jordan to Einstein frame is well defined as long as 
X\ 7^ 2. Equation (37) shows that for x\ ^ 2, if follows that (1 — 2/353) and 



the transformation from the Einstein to Jordan frame is non-singular. For theories 
in which x\ = 2 occurs, this simply represents that the Einstein frame mapping is 
ill-defined and the dynamical analysis should be undertaken in the Jordan frame. 



3.4- Critical points and calculating their stability 

The critical points (also called fixed points) in the Jordan and Einstein frames, x c = 
(xic, x-2c, a^3c) and x c = {x c , y c , z c ) respectively, are the solutions of x[ — x' 2 = x' 3 = 
and dx/d In a = dy/d\na — dz/d\na — 0. These points may be stable, saddle, or 
unstable solutions. In order to study the stability of the critical points we expand 
about these points, x = x c + Ax and x = x c + Ax and consider the eigenvalues of 
the stability matrices M. and M. defined by Ax' = AiAx, and <iAx/rflna = .MAS. 
For any point the matrices M. and M. are related by 



M 



din a 
din a 



TMT- 



d{T) 
din a 



From (40 1, for well-defined transformations (x\ 7^ 2), 



d{T)/d\na — so that we can relate VW, M. and T, 

1 



M = 



1 



(43) 

at the critical points 
(44) 



where T c is the transformation matrix at the critical point. 
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Let the eigenvalues of M at some critical point be Ai, A2 and A3. Then the 
general solution for the evolution of linear perturbations in the JF will be of the form, 

Ax = c ie AlAJV + c 2 e A2AAr + c 3 e A3AAr (45) 

where ci, C2, and C3 are vectors of constants, and N = In a. We require the linear 
perturbations to decay with time so that the critical point is stable. Cosmological 
observations of the cosmic expansion, as measured by the CMB, supernovae and galaxy 
surveys are consistent with a monotonically increasing expansion history such as in 
ACDM [T]-[T7]- We therefore consider stability criteria as the universe expands, taking 
AiV positive, in which case stable fixed points are characterized by negative real parts 
of eigenvalues in the Jordan frame. If the universe had been contracting instead of 
expanding at any point of time, then corresponding fixed points would map to repellers 
instead of attractors. 

If the eigenvalues of M at some critical point (x c ,y c , z c ) are A 1: A 2 and A3, then 
the general solution for the evolution of linear perturbations in the EF will be of the 
form, 

Ax = c ie AlA ^ + c 2 e A2A ^ + c 3 e A3A ^ (46) 

where Ci, £2, and c 3 are vectors of constants and N = In a. We require the linear 
perturbations to decay with time so that the critical point is stable. Our time 
variable in the EF, N, may increase or decrease with the evolution of the universe 
since at a critical point, AN = (1 — 2j3x c )AN and (1 — 2flx c ) could be positive or 
negative. However the eigenvalues in the two frames at a critical point are related by 
A = A/(l — 2f3x c ). The product of these two together hence guarantees that stability 
in one frame implies stability in the other frame irrespective of the sign of (1 — 2(3x c ). 
The sign of the stable eigenvalue in the Einstein frame may therefore be opposite to 
that in the Jordan frame. 

The stability matrix in the Jordan frame, M, is given by: 

Mn = 2/? 2 (3 - 2x ic + AKxi c - 3Kx\ c + 3x 2c + x 3c ) 
3 

(l-2xi c -3Kx 2 lc + x 2c + x 3c ), (47) 



2 

M12 = - 12/3 2 - ^xi c + 6P 2 xi c , (48) 



M13 = - 4/3 2 - ^xi c + 2(3 2 xi c , (49) 
M21 = — + (2 - 4(3 2 )x 2c + 2K(3 - Ap 2 )xi c x 2c , (50) 

TO 

M22 = 3 + 4/3 2 + (2 - 4/3 2 )x lc + "M^lS + (3 _ Ap 2 )Kx\ c 

771 

- 6(1 - 4/3 2 )x 2c - (3-4/? 2 )z 3c , (51) 
= £lc _ mrrxic _ a 

TO TO 

M31 = — + (1 - 4/3 2 )x 3c + 2K(3 - 4(3 2 )xi c x 3c , (53) 

TO 

M 32 - - mrr l Xlc - 3(1 - A(3 2 )x 3c , (54) 

TO 

M 33 = 3 + 4/3 2 + (1 - 4/3 2 )x lc - ^ + + (3 - 4/3 2 )^ 2 c 

TO TO 

-3(l-4/3 2 )x 2c -2(3-4/3 2 ) a; 3 C . (55) 
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where m r = dm/dr, and, 

$F dlnf x 3 z 
f d In <I> X2 y 

The stability matrix in the Einstein frame, M, is given by: 



(56) 



3 

Mn = - 2/3i c + 3x 2 c ~ -i 1 ~ % + Vc + *0, (57) 

M 12 =3f3-^x c , (58) 

M 13 = (3-^x c , (59) 

M21 = — ^ ~ 8f3yc + 6x c y c , (60) 

771 

M22 = - m c - 4pmr f Xc + 3(1 + ~x\ - y c - ~z c ) - 3y c , (61) 

771 z 

~ -i(3x c Af3m r rx c 
M 23 = 1 5 3y c , (62) 

771 ?71 z 

M31 = 6x c z c + 40z e \-l + ~J , (63) 
- 4/3m r r 2 i c 

A4 32 = ^ 3z c , (64) 



M 33 = 4/3£ c [ -1 + - - r c - 3z c + 3(1 + x\ - y c - ~z c ) (65) 

771 / 771^ 



4. Solving the autonomous phase plane equations 

4-1. The critical points 

The seven critical points, PI — P7, are given in Tableland Table[2]in the Einstein and 
Jordan frames respectively. Table [3] gives the Jordan frame fractional matter density 
fl m and effective equation of state parameter w e ft for each of the points, which can 
be compared to observational constraints. 

We will find that r and m, which we can treat as dynamical variables determined 
by the specific theory (the specific form of /(<&)), arc useful parameters to characterize 
the properties and stability of the critical points. 



Note that at the critical points, given that x\ = x' 2 = x' 3 = and (56 1, 



r' = -rxi + lj=0- (66) 

i.e. r is constant at the fixed points. 

At the points P5 and P6 when xi c 7^ and r^O, m and r are related by 

m=-(r + l). (67) 

From ( 28 1 the value of m at these critical points is related to the equation of state by 

m=—^ -. (68) 

3(1 + Weff) 1 ' 

At PI, r = —2 irrespective of the value of (3 or the form of /($), however since 



w e ff = — 1, from (28 1, the coordinates are independent of the value of m. Similarly 



The dynamical viability of scalar-tensor gravity theories 



11 



Table 1. Critical points (PI — P7) of the dynamical attractors for a universal, 
general scalar coupling, j3 in the Einstein frame in terms of dynamical variables 



Point 



y 



pi 

P2 
P3 
P4 
P5 

P6 
P7 



w 

3 



40 
3 



3m 



-4/3+2/3m 



-1 




9-16/3 2 
9 

3/3 2 +m(-9+4/3 2 ) 
4/3 2 (-2+m) 2 



2/3(-l + m) 4/3 2 -8/3 2 m+m 2 (-9+4/3 2 ) 
3m 9m 3 

1 



-(l+m)[-8/3 2 +m(-9+4,3 2 )] 
4/3 2 (-2 + m) 2 



-(l + m)[4/3 2 -8^ 2 m+m 2 (-9+4/3 2 )] 
9m 3 



Table 2. Critical points (PI — P7) of the dynamical attractors for a universal, 
general scalar coupling, /3 in the Jordan frame in terms of dynamical variables 
(x 1 ,x 2 ,x 3 ) 



Point x\ X2 xs 



PI 





-1 


2 


P2 


8/3 2 
4^ 2 -3 








P3 


4/3 
2,3+1 








P4 


16/3 2 


9-16/3 2 





8/3 2 -3 


(3-8/3 2 ) 2 


P5 


3m 


-9m+4/3 2 (-2 + m) 


9m-4/3 2 (-2+m) 


1+m 


16/3 2 (l + m) 2 


16/3 2 (l + m) 


P6 


8/3 2 (l-m) 


-9m 2 +4/3 2 (-l + m) 2 


(l+m)[9m 2 -4/3 2 (-l + m) 2 ] 


3m+4/3 2 (l-m) 


m[3m-4/3 2 (-l + m)] 2 


m[3m-4^ 2 (-l + m)] 2 


P7 


4/3 
2/3-1 
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Table 3. The fractional matter density Q m and effective equation of state w e f f 
in the Jordan frame (JF). For modified gravity theories, observational constraints 
derived from redshift measurements are more directly applicable in the Jordan 
rather than Einstein frame, as discussed in section 12.31 



Point tt m (JF) w eff (JF) 

PI -1 

P2 

P3 I) 



9-4/3 8/3^ 
(3-4/3 2 ) 2 9-12/3 2 



P5 

P6 
P7 



9m 2 +4/3 2 (2-3m + m 2 ) 
8/3 2 (l+m) 2 



3-2/3 
3+6/3 



P4 9-40/3 

U -9+24/3 2 



-8/3 2 +20/3 2 m+m 2 (9-12j3 2 ) 
3m[-4/3 2 +m(-3+4/3 2 )] 

3+2/3 
3-6/3 



at P4, r = and since x% c — 0, from (28 1, again the coordinates are independent of 
the value of m. 

Finally, at points P2, P3 and P7, where x% c — %3c — 0, the coordinates 
{x\, X2, X3} are not explicitly dependent on the specific values of r, m. However, as 



discussed in Sec. 4.3 the stability criteria / eigenvalues at these points are. Note that 



r remains well-defined as given in (56 1 despite X3 and X2, or z and y, asympotically 
tending to zero at a fixed point. 

It should be noted that a specific choice of /($) may only be able to reach a 
subset of fixed points given in this general analysis. For example, consider a scenario 
in which /($) oc <&" (n/ 0), where n is constant, then r = — n at all times, even as 
<f> — > as at P2, P3 and P7. This theory could reach PI only if n = -2, and P4 only 
if n = (for which F($) = and would not give a well defined action). Moreover, 



the consistency relation for P5 and P6 in (67 1 is always satisfied since m = n — 1 for 
this particular choice of /(<&)• 

We briefly summarize the properties of each fixed point below. In each case fl m 
and w e ff are the Jordan frame quantities most easily compared with cosmological 
observations, as discussed in section 12.31 



Point PI: Since £l m = 0, and W e ff = — 1 the point PI can only represent a de-Sitter 
accelerated universe. 

Point P2 : P2 is a generalization of the "</> matter-dominated era" (0MDE) attractor 
in f(R) theories in which the scale factor in the JF has wrong time evolution, a oc t x l 2 
|53j . However for general values of (3, w e ff can be arbitrarily close to zero and a good 
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Figure 1. The effective equation of state uig^ as a function of /3 for the point 
P4. 

approximation to the standard MDE in which a cx i 2 / 3 . P2 could therefore act as 
an approximate CDM era which, for small (3, could be consistent with cosmological 
observations. For \/3\ > v3/2, P2 gives accelerated expansion.. 

Point P3: Since Vl m = at P3, it cannot represent a matter-dominated era (MDE). 
P3 describes an accelerated universe for (3 < —1/2. 

Point P4: For P4 we again have f2 m = 0, so it cannot represent a valid MDE. In Figure 
[I] we show the variation of w e ff with (3 for this point. P4 gives rise to acceleration for 
f3 in the range /3 < — a/3/8, (3 > ^/3/8, or \j3\ < \/3/4. Acceleration has w e // < — 1 
for 13 < -3/4 and (3 > 3/4, and w eff > -1 for \(3\ < For /3 = 1/2, iu e// = 1/3, 

therefore P4 is not accelerated in f(R) theories and general values of (3 thus open up 
avenues for new accelerated critical points. 

Point P5: As m(P5) — > 0, Q, m — > 1 and w e ff —> for all values of (3. Thus 
m(r = —1) = represents a standard MDE. In order for P5 to represent accelerated 
expansion, m(P5) < —1 or to(P5) > 1/2. 

Point P6: Since f2 m = for P6, it can not be used to represent a MDE. P6 gives 
accelerated expansion for a wide range of {m, [3}. In the limit to(P6) — > ±oo the point 
P6 tends to a de-Sitter point with u> e // — > — 1. It is also de-Sitter for m(P6) = 1, 
which is the same as the point PI. As m — ► 1/(1 — 3/4/3 2 ) ornnO", one would get 
acceleration with w e ff — > — oo. 

Point P7: Since f2 m = 0, P7 cannot represent a MDE. For it to represent an 
accelerated universe one requires (3 > 1/2. 

4-2. Comparison with other dynamical analysis of f(R) and coupled CDM theories 

Before discussing the nature of the critical points in detail (which we do in the next 
section) it is useful to make a connection between the solutions that we have found 
and those that have been obtained in the literature, since our general analysis maps 
to specific models that have been considered previously. 
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Table 4. Table showing how theories with purely CDM couplings in which the 
analysis is wholly in the Einstein frame, and f(R) theories with /3 = 1/2, can be 
mapped to subsets of the attractor behavior with a general coupling, /3. 



Point 


CDM coupling, f3 b =0J c = ftl 


f(R) theories, (3= 1/2, [69] 


PI 




Point Pi 


P2 


Point cm 


Point P 2 


P3 


Point d 


Point P 3 


P4 




Point P 4 


P5 


Point 6m (for constant r and to) 


Point P 5 


P6 


Point a (for constant r and to) 


Point P 6 


P7 


Point e 





Solutions for f(R) models with j3 = 1/2 have been obtained in [69] in the JF. 
We can map solutions obtained in [5S] to subsets of solutions found in this paper. A 
summary of the relationship between solutions found here and those in non-minimally 
coupled dark matter and f(R) models is given in Table g 

Our Einstein frame analysis can be applied to CDM coupled dark energy models, 
with fa — 0, Pcdm 7^ 0, e.g. [68] EH 123 [26] . Attractors in these theories have been 

considered with an exponential potential in the EF, V{4>) = Ae^l^ K ^ [6"gll2"4] . 

Only two critical points, which represent a subset of solutions of P5 and P6 with 
Vi oc V2, depend on /1. For these points 

TO = _ ( i + r) = J^. (69) 

We consider in this work a purely scalar-tensor theory of gravity, with no 
higher order metric derivatives, such as the Ricci tensor, R^ v or Riemann tensor, 
RaPiiv Modified gravity theories including such terms have been shown to be 
able to generate cosmic acceleration, for example Gauss-Bonnet gravity /(G) where 
G = R 2 — ARVR^v + R a @v 'Rap^v (which is motivated by effective low-energy actions 
in string theory) H2 E2J E01 1^1 EH EH E51 ES] ■ Our analysis is not applicable 

to higher derivative gravity theories. Dynamical analyses [62j [66] have been used to 
show, however, that f(G) models are highly constrained by cosmological data. 

4-3. Stability criteria for the dynamical attractors 

We present here the criteria for stability in each of the critical points discussed in 
Sec. AA_ and shown in Tables [l] and [2] A viable cosmology requires an extended 



matter-dominated era (MDE) followed by late time acceleration in the Jordan frame. 
Thcrfore the MDE needs to be a saddle point and the acceleration can be stable or 
saddle. The stability criteria presented here is for an expanding universe. In the case 
of contraction, stable attractors map to unstable repellers and vice versa. 



In Sec. 4.1 we discussed the dependency of the observables f2 m and w e ff on to and 
j3. At certain critical points, PI, P4 — P6, r is a constant or equal to — (m+ 1), while 
for the points, P2, P3 and P7, r is unconstrained. Here we find that the conditions 
for saddle matter domination and stable acceleration can be expressed in terms of 
conditions on to, r and the coupling, (3 for all points PI — P7. 
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m and r are considered variables determined by the specific theory, in terms of 
which the stability of the critical points is described. We denote the value of m at the 
critical point Pi by m, = m(JPi) and the derivative of mi with respect to r at that 
point by mi r = dm/dr(Wi). In this section we consider conditions to give saddle MDE 
and stable acceleration in the Jordan frame. We then give conditions on the specific 



evolutionary paths for m and r that yield viable cosmologies in Sec. 4.4 and discuss 
conditions for saddle acceleration in Sec. 14.51 

Point PI: 



The eigenvalues for PI are, 




where m\ = m(r = —2). For the point to be stable we need the real parts of all 3 
eigenvalues to be negative (or zero). This amounts to < m(r = —2) < 1 irrespective 
of the value of (3. PI is a saddle point otherwise. 



Point P2: 

The eigenvalues for P2 are, 

1 3 1 

~ 2 ~ 3 - 4/3 2 ' m 2 (3-4/3 2 ) 



9to 2 - 4/3 2 {m 2 (-l + 2m 2 ) + m 2r r(l + r)} 



±Af3 2 \Jm\ + 2m\ + m i I 1 + 2m 2r r(-l + r)} - 2m 2 m2 r r(l + r) + m 2 2r r 2 {l + r) 2 

As an illustrative example, if m 2 is constant so that m 2r = 0, the eigenvalues reduce 
to, 

_1_ 3 9m 2 +4/3 2 (2-m 2 ) 
2 3-4/3 2 ' ' m 2 (3-4/3 2 ) 

Since at least one eigenvalue is positive, P2 can either be a saddle point or unstable 
and does not represent a stable acceleration point even though w e // < ~ 1/3 for 
\(3\>VS/2. 

Point P3: 



The eigenvalues for P3 are, 
3 + 2/3 2 

IW m 2 (l + 2/3) [~^ + ^ r(1 + r) + 3m ' (1 + ^ 
±P^Jmj + 2ml + + 2m 3r r(-l + r)} - 2m 3 m 3r r(l + r) + m\ T r' 1 {\ + rf 

For example, if m 3 is constant so that m 3r — 0, the eigenvalues reduce to, 
3 + 2/3 6 + 8/3 6m 3 +4/3(-l + TO 3 ) 
1 + 2/3' 1 + 2/3' m 3 (l + 2/3) ' 

In the range —3/4 < (3 < —1/2 the first and second eigenvalues are negative, in this 
constant m 3 case. To get stable acceleration, one must find corresponding values of 
m 3 so as to get the third eigenvalue negative or zero as well. 



Point P4: 
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The eigenvalues for P4 are, 

o 3 16/? 2 (l + m 4 ) 



8/3 2 ' m 4 (3-8/3 2 ) ' (?0) 



where 7724 = m(r = 0). Accelerative expansion occurs for (3 < — -y/3/8, f3 > -\/3/8, or 
|/3| < \/3/4. The first eigenvalue is negative, and the second one is negative or equal 
to zero for (3 < —3/4, (3 > 3/4, or \(3\ < ^/3/8. In the overlapping regions, i.e. when 
/3 < —3/4, (3 > 3/4, \(3\ < \/3/4 we have almost all conditions for stable acceleration 
satisfied. We can then find corresponding values of 777,4 so that the third eigenvalue 
also becomes negative or zero. For \(3\ < a/3/4 and -1 < m(r = 0) < 0, P4 always 
represents a stable accelerated era with w e ff > — 1, with no condition on m^ r . The 
point P4 is therefore of interest and we discuss it further in sections |4.4| and |4.5| 

Point P5: 




The general expression for the eigenvalues of P5 is non-trivial. The eigenvalues of P5 
in the limit of a matter dominated era (with \m^\ << 1) are approximately, 

3(1 + m 5r ), 

Models with 7775 < are not acceptable since the eigenvalues diverge as 777,5 — > 0~, 
so the system would not remain at P5 for sufficient time to give rise to a viable CDM 
dominated era [53] ■ In order to get a valid saddle MDE as 777,5 — > + we need, 

m(r < -1) > 0, m r (r < -1) > -1, m(r = -1) = 0. 

We find that P5 is never stable in cither of the two regions that give accelerated 
expansion (777,5 < — 1 and 7775 > 1/2), for values of [3 approximately in the range 
|/3| < 0.86. For (3 outside of this range P5 can give rise to accelerated expansion, 
however these values of (3 are generally disfavored by observations (24]. Hence P5 
best describes a MDE. 

Point P6: 



The eigenvalues for P6 are, 

-Af3 2 + 8/3 2 m 6 + m 2 (9 - Af3 2 ) -8/3 2 + 12/3 2 m 6 + tt7 2 (9 - 4/3 2 ) 
m 6 [-4/3 2 + m 6 (-3 + 4/3 2 )] ' T77 6 [-4/? 2 + m 6 (-3 + 4/3 2 )] ' 

8/3 2 (-l+ 777 2 )(l+ 777 6r ) 

777 6 [-4/3 2 + m 6 (-3 + 4/3 2 )]' 1 ' 

We note that there is a symmetry with respect to positive and negative values of (3 
since w e ff and the eigenvalues are all functions of f3 2 . 

In order to characterize and study the properties of P6 in more detail, we limit 
our analysis to (3 in the range \(3\ < 1/2, in which we find that it is accelerated 
and stable in five distinct ranges. This choice of (3 is well-motivated, the value of 
(3 is constrained by CMB observations with \(3\ < 0.1 in [23] and encompasses f(R) 
theories with (3 — 1/2. Hence |/3| < 1/2 represents a broad and interesting range for 
viable scalar-tensor theories. 

There are five regions of stable acceleration in P6, labelled (A)-(E), which depend 
on the magnitude of m at the fixed point, 77i6, and the sign of the derivative dm/dr 
at the fixed point, written m§ r . 

[I] 777 6r > -1 
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(A) m 6 < min [-1, 2|/3|/(2|/3| - Vs)}: w eff > -1, 

(B) 4/3 2 /(4/3 2 - 3) < m e < 0: w eff < -1, 

(C) m 6 > 1: w e// < -1, 

[II] m 6r < — 1 

(D) 2|/3|/(2|/3| + V3) < m 6 < 1: w eff > -1 

(E) -1 < m 6 < 2|/3|/(2|/?| - y^): This point exists only for (3 in the range 
\(3\ < n/3/4 for which w e // > -1. 

Region (E) does not exist for f(R) theories, so as with P4, general values of (3 
therefore open up avenues for new accelerative critical points. We discuss this region 
further in sections 14.41 and 14.51 

Point P7: 



The eigenvalues for P7 are, 
3 — 2/3 2 

1^9' m 2 (l - 2/3) P 7 ~ + T) ~ 3m ' ( ' 1 + P) 

±(3\J ™| + 2m 7 + to 2 {1 + 2m-[ r r{— 1 + r)} — 2m 7 m 7r 7 , (l + r) + rnj r r 2 (l + r) 2 

For example, if m 7 is constant so that m 7r — 0, the eigenvalues reduce to, 

3 - 2/3 6 - 8/3 6m 7 + 4/3(1 - m 7 ) 
1-2/3' 1-2/3' m 7 (l - 2/3) 

Stable acceleration is able to occur for the range 1/2 < (3 < 3/4, where the first and 
second eigenvalues are negative (in the constant m-; case), and a suitable value of m 7 
can make the third eigenvalue to also be negative or zero. 



4-4- Classification of F(<fr)R scalar-tensor theories 

We classify general F(Q)R scalar-tensor theories, with values of (3 in the range 
\P\ < 1/2, on the basis of trajectories of the m(r) line on the (r, m) plane since such a 
trajectory completely specifies a theory's dynamical evolution as the universe expands. 
These trajectories depend on the functional form of F(Q), or from an alternative 
perspective on the value of [3. We therefore consider conditions on viable cosmologies 
as conditions on these trajectories. 

• We assume that the function / is a C°° function (i.e. it is differentiable 
for all degrees of differentiation). Also / and all of its derivatives are non- 
singular and single valued. For the conformal transformation and scalar field 
redefinition to be valid we require that F($) > 0. This implies that f(Q) is 
monotonically increasing and single valued for all As such r(<I>) and to(<I>) 
are also single valued. Therefore m(r) is single valued and we can consider 
evolutionary trajectories in terms of the function m(r). 

• As mentioned earlier, a viable cosmology requires an extended matter dominated 
era followed by late time acceleration. The matter dominated era therefore needs 
to be a saddle point. From our detailed analysis in Sec. |4.3| we find that only 
the point P5 with 7715(7' = —1) — > + can be used as a standard saddle matter 
dominated point with J7 m = 1, w e // = and a cx t 2 ^. This point is denoted as 
P5W . 
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• The accelerated expansion needs to be a stable or saddle point. Valid acceleration 
points (for < 1/2) are PI, P4, and P6. Transitions from P5 (0) to PI and 
regions (A)-(D) of P6 are qualitatively similar to those in f(R) theories [55] . 
The point P4 and the region (E) of point P6, however, are not present in f(R) 
models, we therefore discuss their properties in detail below. Interestingly, both 
of these exist as stable accelerated epochs only for |/3| < -s/3/4 and have very 
similar properties. 

• The points P5 and P6 lie on the critical line, to = — r — 1. Valid m(r) trajectories 
should not enter a critical point from a forbidden direction. For the point P5^ ' 
to exist we require m 5r > — 1. We can therefore mark forbidden direction regions 
on the critical line around the point PS^ -*. Forbidden direction regions also exist 
around P6, with transitions to the regions (A), (B) and (C) of point P6 requiring 
that m§ r > — 1, and m§ r < — 1 for transitions to regions (D) and (E). The point 
P4, by contrast, has no such forbidden direction regions around it. 

The characteristics of the 5 general classes of scalar-tensor theory trajectories that 
give stable acceleration are summarized below. 

I: The m(r) curve does not connect the standard matter dominated point to the 
accelerated attractor solutions. Models of this class either bypass the matter era 
and directly fall onto an accelerated attractor or go through a </>MDE. For general 
values of (3, however, w e ff for the point P2 can be arbitrarily close to zero and 
give a standard MDE instead of a </MDE. 
II: The m(r) curve connects the standard matter era P5^ ' to PI or asymptotically 
(r — > ±oo ) to P6, giving a stable accelerated de-Sitter expansion. 

Ill: The m(r) curve intersects the critical line in the region (B) of point P6. These 
models have very short matter dominated phases which would not allow sufficient 
time for structure formation. The resulting acceleration after matter domination 
has w e ff <C — 1. 

IV: The m(r) curve connects the standard matter era P5(°) to the region (D) of P6 

which has acceleration with w e ff > — 1. 
V: Models of this class connect the standard MDE P5(°) to the point P4 in 
— 1 < m(r — 0) < 0, or to region (E) of P6. These models exist only for 
\j3\ < Vo/4, hence P4 and region (E) of P6 are not accelerated for f(R) theories. 
The final acceleration has w e ff > — 1. 

Classes I-IV are qualitatively similar to those in f(R) theories. Class V is not 
open to f(R) theories and we discuss this class here in more detail. 

Class V covers the transition between points P5' - 1 and P4 or region (E) of P6. 



As shown in Sec. 4.3 the point P4 and region (E) of P6 represent accelerated epochs 
with w e ff > —1 for |/3| < a/3/4. P4 is stable for —1 < m(r = 0) < irrespective 
of the value of m± r , and region (E) of P6 (which has —1 < mg < 2|/3|/(2|/3| — V%)) 
is stable for mg r < — 1. If we can find a function /(<&) which takes us from to 
cither of these points then this would represent a valid trajectory for cosmic evolution. 

In Figure [2] we show possible trajectories on the (r, to) plane for a Class V model. 
Consider the case when Pb^ is connected to P4 by a straight line. We require the 
slope of this line to be between and -1 so that it can intersect the TO-axis between 
and -1. The resulting equation for the m(r) curve that we are looking for is, 

m=--(l + r) (72) 
n 
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Figure 2. Class V models of scalar-tensor theories on the (r, m) plane. The 
solid line with a slope of -1 is the critical line m = — r — 1 and the dashed curves 
are possible trajectories for a Class V model. The point (—1,0) is the 

starting saddle matter dominated point, with the triangles around it representing 
the forbidden direction regions. The m(r) curve for a Class V model intersects 
either the m-axis between and -1 corresponding to the point P4, or the critical 
line in region (E) of point P6, to get a stable accelerated epoch. 



where n > 1. Combining (72) with (17 1, (28 1 and (56), the resulting second-order 
differential equation in /(<&) gives solutions for Class V models of the form, 

/($) =a($^ +b)^ (73) 
where a and b are constants, and a > 0. This gives, 

(74) 



In order to have the correct sequence of matter domination and acceleration we need 
r to be increasing with time, i.e. dr/da > 0. Using $ < and a > this corresponds 
to b > (here dots denote derivatives with physical time in the JF). In f(R) theories, 
<f> < represents a curvature, R 7 that decreases with time consistent with observations. 

The above form of /(<&) can be used to calculate the mass of the scalar field $ in 
these scalar-tensor theories [35] , 

ml = = -f- f — + $ - 4/($)J (75) 

For large values of 4>, b > corresponds to to^ < 0. Thus this simple /(4>) leads 
to scalar fields with imaginary mass. If we have some general curve joining P5<W and 
P4 we can break it up into an infinite number of straight lines with at least some being 
of the form of the above /(<&) (with slope between -1 and 0) and again argue that it 
would lead to an imaginary mass for the scalar field. A transition from to region 

(E) of P6 will similarly need an imaginary mass scalar field because this transition 
also goes in the direction of increasing r with a negative value of m. Hence Class 
V have a matter dominated era followed by an accelerated epoch with w e tf > — 1, 
but relate to theories with imaginary scalar field masses in the Einstein frame. In 
the specfic case of f(R) models this translates to imaginary mass scalar fields arising 
when f RR < - see Refs. [5SJ [55] for details. 
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4-5. Criteria for viable models with saddle acceleration 

Till now we have only studied viable cosmologies which go from a saddle matter 
dominated era to a stable accelerated universe. Even though we know that dark 
energy is the dominant component today we can not be sure of what the future holds. 
Therefore it may be possible that acceleration is not a permanent feature of our 
universe [87] . 

We first consider acceleration that is stable on the subspace z = ry, on which 
saddle point CDM attractors are stable, where r (and hence m) is constant. The 
evolution equations written earlier reduce to, 

x 1 = [3[1 - x 2 - (1 + r)y] + 2(3(2 + r)y - - x 2 + (1 + r)y] (76) 

y> = - A(3rxy { — ) - 8f3xy + 3y[l + x 2 - (1 + r)y] (77) 
\m J 

where r = —(1 + m). The solutions of these evolution equations are the points P2, 
P3, P5, P6, and P7. Since m is now constant, the points PI and P4 are just special 
cases of P6 with m = ±1. 

For |/3| < 1/2 we see that the only points that can be accelerated are P5 and P6. 
It can be shown that P5 is not stable (even on the subspace z = ry) in the regions 
that it is accelerated. The point P6 ( m ~^ m \ 4/3 2 -8/3 2 m W (-9+4/3^ has Q m = , 
an effective equation of state, 

-8/? 2 + 20/3 2 m + m 2 (9 - 12/3 2 ) 



3m[-4/3 2 + m(-3 + 4/? 2 )] 



(78) 



and its eigenvalues are given by the first two in (71 1, 



-4/3 2 + 8/3 2 m 6 + m 2 (9 - 4/3 2 ) -8/3 2 + 12/3 2 m 6 + m 2 (9 - 4/3 2 ) 

(79) 



m 6 [-4/3 2 + m 6 (-3 + 4/3 2 )] ' m 6 [-4/3 2 + m 6 (-3 + 4/3 2 )] 

P6 is accelerated (w e ff < —1/3) and stable on the subspace z = ry in the following 
three regions: 

(A) m 6 < 2|/3|/(2|/3| - ^3): w eff > -1 

(B) 4/?7(4/3 2 - 3) < m 6 < 0: io e // < -1 

(C) m 6 > 2\/3\/(2\P\ + V3): w e // > -1 for m 6 < 1, and w eff < -1 for m 6 > 1 

Notice that there is no condition on m,Q r , therefore there is no forbidden direction 
region around P6. The first graph in Figure [3] shows a possible trajectory of the m(r) 
curve, connecting P5 1 - ' to either of the regions (A) and (C). We can also connect 
to the region (B) of P6. It is likely that all scalar-tensor theories which connect to 
the regions (A) and (B) are tachyonic because of reasons discussed earlier (since they 
have negative values of m and are moving in the direction of increasing r) . This only 
leaves region (C) to be a viable cosmology. Theories which connect P5(°) to region 
(C) of P6 are therefore possible models for the evolution of the universe from a saddle 
matter dominated era to a saddle accelerated expansion. 

On a general subspace (away from z — ry) the point P4 admits acceleration for 
|/3| < a/3/4, not compatible with f(R) theories, that gives rise to a second saddle 
acceleration attractor. If we consider positive values of m(r — 0) then we see that 



the third eigenvalue in (70 1 is positive but the first two are still negative. Thus P4 
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Figure 3. Scalar-tensor theories which lead to a saddle accelerated expansion of 
the universe. The solid line with a slope of -1 is the critical line m = —r — 1 and 
the dashed curves are possible trajectories for a scalar-tensor theory. The point 
P5<0) (—1^0) is the starting saddle matter dominated point, with the triangles 
around it representing the forbidden direction regions. In the first graph the m(r) 
curve connects the matter dominated point to one of the three accelerated regions 
(A, B or C) of the point P6. The accelerated expansion is stable on the subspacc 
z = ry. In the second graph the m(r) curve connects the matter dominated point 
to the point P4 on m(r = 0) > which gives a period of, potentially transient, 
acceleration. 

with m(r — 0) > represents saddle acceleration on such a subspace, for \/3\ < v3/4. 
The second graph in Figure [3] shows a possible trajectory from the matter dominated 
point P5'- -' to the saddle acceleration point P4. 

5. Conclusions 

In this paper we have established the dynamical attractor behavior in scalar-tensor 
theories of dark energy, presenting a complete, consistent picture of evolution in both 
the Einstein and Jordan frames. We discuss critical points for the dynamical evolution 
and show how in the two frames the stability matrices for these points are related by 
a similarity transformation. 

We carry out a general analysis for values of the coupling (3 in \j3\ < 1/2, and 
find that there are five classes of evolutionary behavior, of which four classes are 
qualitatively similar to those for f(R) gravity. The fifth class only exists for values of 
(3 in the range |/?| < \/3/4, i.e. it is not present for f(R) gravity. This class of models 
has a standard matter dominated phase followed by acceleration with w e ff > — 1, but 
for stable acceleration seems to always have tachyonic scalar fields. 

We then relax the condition of a stable accelerated universe and study transitions 
from the allowed MDE to a saddle accelerated universe, again for < 1/2. We 
find that on the subspace z = ry (r = constant), which contains the suitable MDE, 
the point P6 is accelerated and gives evolution without tachyonic scalar fields if 
m6 > 2|/3|/(2|/3| + \/3)- This class of behavior is qualitatively true also for f(R) 
theories ((3 — 1/2). In addition there is a possibility of going from the MDE to the 
saddle non-phantom acceleration point P4, only open to \(3\ < \/3/4. Therefore if we 
allow the current accelerated expansion of the universe to be temporary then there 
are a much broader variety of f(R) and F(<fr)R theories that are dynamically valid. 
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Our results open up interesting questions about broader dark energy theories. In 
particular, in this paper we have used a particular form of the potential in the EF 



given by Eq. (18 1 appropriate for a wide class of theories, including f(R) theories, in 
which the Jordan frame action has no scalar potential. It will be interesting to see 
the implications of using a more general form of the potential. 
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